Pion Correlations and Two-Particle Bound States 



M.Shmatikov 

Institute for Theoretical Physics II, 
Ruhr Universitdt Bochum, Bochum, Germany 

2£ i and 

Q\ [ Russian Research Center "Kurchatov Institute", 

123182 Moscow, Russia f] 

< 



=5 



Abstract 



; A method to investigate possible loosely bound dibaryon states is suggested. 

The momentum dependence of a two-pion correlation function is shown to 
depend on the (non) existence of a bound state. 

\D ' 

in 
o 

Tj" | The only dibaryon state which existence does not cause at present any 

doubts is the deuteron. There is also a plausible candidate suggested for the 
explanation of the forward angle cross section of the pionic double charge 
exchange [|lj. It is unnecessary to stress that its experimental observation 
became possible because of the extremely small width of the resonance: 
T < 0.5 MeV g. 

Other dibaryon states with much larger width, exemplified by a AA 
bound state 0, are also conceivable. A naive estimate of its width r > 2 Ta ~ 
300 MeV prompts, however, that it will hardly manifest itself as a resonance 
| peak. 

Here we address the problem of observing possible bound states of two- 
baryon resonances. To be more specific we consider a pair of two baryon 
resonances (Bi and B 2 ), their interaction being strong enough to produce a 
(loosely-) bound state. Were the bound-state components stable it would, 
like the deuteron, exist for an infinitely long time. However, because of the 
baryons' decay the bound state itself exists for a finite time. The only signal 
of a bound state emergence will be then the localization of the B\ and B 2 
particles in some finite spatial region. 

A dominant decay channel of a baryon resonance is usually that with the 
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7r-meson emission: 

B -> B' + vr 

Then the (Si S2) bound state will decay with the emission of two pions: 

{B 1 B 2 ) -> + B' 2 + vr + Ti (1) 

In the case of a (AA) bound state the decay pattern ([[]) translates into 

(AA) — >iV + iV + 7r + 7r (2) 

Hereafter we shall address specifically the process (0) though all the results 
are applicable to a more general case dH). 

Spatial distribution of two A-isobars can be scanned using the two-particle 
correlation in the momentum space. The latter in the form of the pionic in- 
terferometry is routinely applied for extracting information about the space- 
time structure of the emitting sources in heavy-ion collisions ||. 

Forbidding complexity of the multiparticle process ([!]) (and (|2|)) necessi- 
tates making some simplifying approximations. In the vein of the N c — ► 00 
approach [[| we consider the A-isobar (and the nucleon) as a heavy source of 
7r-mesons. Then the process (Q) can be viewed as an emission of the pions by 
two heavy sources superimposed by the motion of the sources which is not 
affected by the pion emission. 

Separating the system under consideration into two components: a fast 
(A — > Ntt decay) and a slow (relative motion of two A-isobars and nucleons) 
ones we apply the adiabatic approximation 

Consider first two A-isobars at a distance -Raa- Correlation between two 
(identical) pions due to the quantum mechanical uncertainty depends on the 
sources' separation -Raa- In what follows we use the formalism developed in 
|J. Let ip n7T be the wave function of the tttt system. Its time evolution is 
governed by the hamiltonian H: 

4>7m(r a , r b ; t) = exp (-iH(f a , f%)(t - t )) Vw(r a , r b ; t ) (3) 

We define also an asymptotic state with definite values of the 7r-meson mo- 
menta 

4>p a ,p b {r a , n; t) = exp (-i(E a + E b )t) exp (i(p a f a + p b r b )) (4) 
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where E a>b = Jm 2 + p 2 b are the particle energies. Then the measured two- 
particle amplitude equals 0: 

A(Pa,Pb) = Hm^oo J d 3 r a d 3 r b ip^^r a , r b , t ) 

x [exp [-iH(f a , r b ){t - to)] fe aA (f , r b ; t)]* 

The values of t for which this expression is valid will be discussed below. 

It is convenient to introduce center-of-mass and relative coordinates of 
two 7r-mesons: 

r = f a - f b ; R = Ta - n (6) 

and also the total and the relative momenta of the pair: 

q = Pa~Pb; P = Pa+Pb (7) 

Using these coordinates the asymptotic wave function can be rewritten as: 



'Pa,Pb 



(f a ,f b ;t)=e^ PR - Et )e^ (8) 



where E = E a + E b is the total energy of the tttt pair. The total energy E 
may be separated then into the c.m.s motion energy E cm and the relative 



motion energy E Te \. The former obviously reads E cm = y P 2 + Am 2 . Then 
E Te \ = E — E cm . In a frame where P ~ the relative motion energy reduces 
to E Te i ~ q 2 /(Am v ). We make a simplifying assumption that the emitted tt- 
mesons do not interact with each other. (The validity of this approximation 
is discussed below.) Then the two-pion amplitude (^) can be cast in the form 



A(p a ,p b ) = e iEt °J d\SRe- iM e^ 2 ^ (r+^R- 7 --^) (9) 

We assume that two pions are emitted in the decays of baryons indepen- 
dently, implying that at some time moment t a the two-particle wave function 
ipTm factorizes: 

ip7m{r a , n; t a ) = -j= [ip a {r a ] t a ) ^&(f 6 ; t a ) + ^ a (f 6 ; t a ) tp b (f a ; t a )} (10) 

Consider the case when the 7r-meson, in the state described by the ip b 
wave function, was emitted earlier than that in the ip a state at some time 
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t b < t a . Then by the moment t a (when the second pion was emitted) its wave 
function has acquired some phase: 

Mr; t a ) = e-^f-fc) Mr; U) (11) 

where h is the one-particle Hamiltonian (of free motion): h(r) = m — ^V^. 
The two-pion wave function reads then ||: 

M(r a , n; t) = [ Mr a ; t a ) e -<^)(*»-*») Mn; h) 



+Mn; Q e- ih ^^Mf a] t b ) 
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The time evolution of the wave function(s) is described the most easily by 
introducing the Fourier decomposition into the momentum eigenstates. The 
two-particle amplitude corresponding to the appearance of the second n- 
meson at the time moment t a can be represented then as follows ||: 

A$(p a ,p b ) = V2e(ta-t b )J0,0,d"rd^Re^- p y n e^^ 

(IS) 



x cos (f ) % (f + k; t a ) i> b (I - k; t b 



Here ip(k; t) is the Fourier transform of the one-particle wave function defined 
as follows 

mt)= J ^-J(k;t)e^~^ (14) 



and to = ym 2 + k 2 . In eq.([L3|) Q = k a + k b and k — (k a — k b )/2. Finally, 
f2(Q, k) = E — to (Q/2 + k) — oj [Q/2 — kj. The integral over the c.m.s. 

coordinate R yields (2%) 3 8^ (Q — P), while the integration over r produces 
a Fourier image of the cos(<fr*/2) function $(<f/2, k): 

Aa a b \p a ,p b )= e(t a -t b )f^^ p / 2 ^$(q/2,k) 
xMP/2 + k;t a )MP/2-k;t b ) 

An alternative situation when the particle in the state a was emitted earlier 
at the time t a and the other particle appeared later at the time t b yields 
similar expression M: 

A a b b \pa,p b )= 9(t b -t a )f^e^ p ^$(q/2,k) 
xMP/2 + k;t a )MP/2-k;t b ) 
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The total amplitude is the sum of two components (eq.([H|) and flT6|)) corre- 
sponding to two time orderings of the particles emission: 

A ab (Pa, Pb) = A$ (pa, Pb) + (fa, Pb) (17) 

Using this amplitude we calculate the two-particle cross section which is 
defined as follows || 

P2(Pa,Pb)= J dt a dt b dt a > dt h > p aho!h > A* a , h ,{p a ,p h )A ab {p a ,p b ) (18) 

ab,a'b' 

Here p ab a iy is the density matrix of the source which contains information 
on the probability distribution for the two-particle quantum numbers (a, b) 
and for the emission times (t a ,t b ). 

Independent emission of two particles is ensured by a factorizable ansatz 
for the density matrix ||: 

Pab,a'b' = U aa' P(*a, V) P bb' Pfa, fy) (19) 

The probability of the two-particle emission consists of four terms: 

P 2 (p a ,p b ) = P^ + + + P^ (20) 

The reader is referred to § for the proof that P (ab) = P (kl) = and P (aa) = 
P^ so that 

P 2 (p a ,Pb) = 2P {aa) (21) 
Following || we define the single particle density S(X, K) of the source: 

S{X, K) = J d 4 x e^ K ^ p{X° + x°/2, X° - x°/2) 

x E aa > v ao > MX + x/2) r n > (X - x/2) W 



where 



2 

and the 4-momentum K equals: 



X° = \(*. + tj), *" = «„-«„• (23) 



j? = A±ft; A„ = ^ (24) 
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Then the two-particle probability (|T8| ) can be expressed in terms of the single- 
particle density 0: 



P2(p a ,Pb) = fd 4 Xd 4 Y [S{X,K + q/2)S{Y,K-q/2) 
+e *<i(x-Y) S (x 7 K)S(Y,K) 



(25) 



Here q = E a — E b and q = p a — Pb- Remind that the expression for the two- 
particle probability was obtained in the frame where p a + p b ~ implying 
that go ~ 0. Single particle spectrum is given by the expression: 

P 1 (p)= Jd 4 XS(X,p) (26) 



with the 4-momentum p = (yp 2 + m%,p) 

We are now in position to investigate the two-particle correlation function 
C(p a ,p b ) = C(q,K): 

c ^~-im§» < 27 » 

We simplify this expression by separating the time and the spatial de- 
pendencies of the single-particle density S(X,K) (eq.(|22|)). To this end we 
introduce a short-hand notation 

p(T,t) = p(X° + x°/2;X° -x°/2) (28) 

denoting the rest of the S(X, K) as S: 

S(R, K) = [ d 3 r e iJ?f £ "aa'A (R + r/2) (R - f/2) (29) 

aa 

so that 

S(X, K) = Jdt e lEt p(T, t) S(R } K); X= (T, R) (30) 

Then the two-particle probability (pop can be cast in the form: 



P 2 = J dT x dT 2 d 3 R 1 d 3 R 2 J dt x dt 2 {e iE ^ e lE ^ p x p 2 S{R x ,p a ) S(R 2 ,p b ) 

+e i(E a -E b )( Tl -T 2 ) e i^R : -R 2 ) e iK 0tl e iK t 2f)i ^ g^ g^ j 

(31) 
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where we introduced for notational brevity pi = p(T i; ti); i = 1, 2. Note that 
in the considered frame K w the energies of the pions are almost equal to 
each other (E a « E^), implying that K « -E a (~ -£?&)• The P2 probability 
may be simplified further by introducing an (energy-dependent) factor 

X(E) = J dTdtp{T,t)e lEt (32) 

We arrive at 

P2(p a ,Pb) = X(E a ) X(E b ) J d 3 Ri d 3 R 2 

The single-particle probability in these notations reads: 

P 1 (p a ) = X(E a ) Jd 3 RS(R,p a ) (34) 

with the similar expression for Pi($,). 

We obtain finally the two-particle correlation function: 

C(Pa,Pb) = 1 + 

/ d 3 R 1 d 3 R 2 e (A-A)(5i-A) 3 (R x , § fa, (35) 

{fdZRxSiR^Pa) Jd 3 R 2 S(R 2 ,p b )y 1 

Let us consider now a bound state of two baryons. As it was stated 
above, in the N c — > 00 limit they can be considered as heavy sources of tt- 
mesons. Neglecting recoil due to the pion emission we can write the spatial 
distribution of the sources as 

^ = 7 {5(r - R AA /2) + 5(f+ Raa/2)} (36) 

where -Raa is the distance between two A-isobars (baryons) and 7 is a nor- 
malization factor. (Hereafter for notational brevity we omit the wave function 
subscript denoting the charge state of the emitted pion). Inserting the wave 
function into the expression for the S (eq. (p9|) ) we arrive at 

S(R,p) = 7 2 {5(R-R AA /2) +5(R + R AA /2)} (37) 
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and the corresponding single-particle distribution 

P 1 (p a )=X(E a )2 1 2 



(38) 



The two-particle distribution calculated with the same source function (|3"7j) 

reads 

P2(p a ,Pb) = 27 4 (1 + cos(^Raa)) (39) 

where q = p a — pt,. Note that in deducing P 2 only those configuration were 
taken into consideration when each heavy baryon emits a pion, i.e. the tt- 
mesons are emitted in different spatial points. 

Inserting the single- (eq.([I3|)) and two-particle (eq. (^)) distributions in 
the expression ([27]) we arrive at the correlation function for pions emitted by 
two heavy sources separated by the distance R (since it will not cause any 
confusion the subscript A A is omitted hereafter): 

C n = 1 + cos(qR) (40) 

Again, only those configurations were taken into account when each A-isobar 
de-excitates by the emission of the 7r-meson. It implies that the correla- 
tion function is contributed by the configurations when the emission points 
of the 7r-mesons are spatially separated. Note also that since p a +Pb ~ 0, the 
relative momentum of two pions equals q ~ 2p a . 

The correlation function (|40"D depends parametrically on the distance be- 
tween pion-emitting sources R. Two baryons being bound, their spatial dis- 
tribution is governed by the bound-state wave function <ft(R): the probability 
that the baryons are at the distance R is given by the \<fi(R)\ 2 . The spatial 
distribution of pion-emitting sources does not affect the single-particle dis- 
tribution (P8|), whereas the two-particle distribution ( |3"§D and, consequently, 
the two-pion correlation function is modified. The observable correlation 
function 

C = J d 3 RC n \(j)(R)\ 2 (41) 

reads 

C(p a ,p b ) = 1 + Jd 3 R cos(qR) |0(i?)| 2 (42) 

This expression shows that the two-pion correlation function with q varying 
"scans" the wave function of the baryon relative motion. 
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Eq. (O) becomes more "transparent" in the case of an S'-wave two-baryon 
bound state. Assuming a simplest deuteron-like wave function 

«*> = Vshr (43) 



where k is the bound-state momentum (k = ye M, e is the binding energy 
and M is the reduced mass of the pair of baryons) after straightforward 
integration over angular variables in eq.(PE3) we arrive at 



C = l + 2. (e-^ ^^dR (44) 
J qR 

Introducing the \ = ?/2k ratio and calculating the integral we obtain the 
two-pion correlation function reading: 

c = i+ arctan X 
X 

The dependence of the correlation function C on x is plotted in the figure. 
Let us consider it in more detail. Given a pair of bound A-isobars they are 
localized in a domain with a spatial size R ~ l/n. For the values of q small 
compared to 1/R (the ratio % — <?/2k <C 1) both sources emit the 7r-mesons 
coherently resulting in the 100%-interference and, hence, in C — > 2. In the 
opposite case (g > 1/E) because of the oscillation of the relative phase (see 



eq. (|4"2|) ) the interference between the emitted pions is smeared out. It results 
in C — > 1 at large values of x- The devitation of the correlation function C 
at large x' s from the unity is proportional to the bound state momentum: 

o-*~l-7 (46) 

In the case of unbound A-isobars the wave function of their relative mo- 
tion (squared) does not have the sense of the probability of finding two par- 
ticles at some distance R. (It is reflected, in particular, in the different 
normalization of the wave function of the continuous spectrum). This formal 
obstacle can be circumvented by considering a continuous-spectrum wave 
function as that of a bound state with vanishing binding energy. Then we 
get from the eq.fllSD that C = 1 for all values of q. This conclusion fol- 
lows immediately from the qualitative consideration. Indeed, since in the 
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Figure 1: Momentum dependence of the correlation function in the presence 
(solid line) and absence (broken line) of a bound state 

continuous spectrum arbitrarily large distances between emitting sources are 
possible the relative phase controlled by the qR product will be large for any 
(nonvanishing) value of q. Inspection of (f42"D shows that because of rapid 
oscillations of the integrand C = 1 for any unbound state. 

When deducing correlation function ( fl2|) several approximations have 
been made. The first approximation is that the emission of a 7r-meson changes 
the discrete numbers of a baryon B only (e.g. A — > N) and does not affect 
its position. Such an assumption seems to be justified in the N c — > oo limit 
where baryons may be considered as (infinitely) heavy particles. Before pro- 
ceeding to the discussion of other approximations made we consider in more 
detail (possible) bound states. The very concept of a bound state of two 
baryons B\ and B 2 assumes that they do not lose their identity. It can be 
translated in the requirement that the characteristic distance R between two 
particles exceeds their doubled size b: R > 2 b. Taking as a benchmark of 
the baryon's size that of the nucleon {b ~ 0.5 fm) we arrive at the con- 
clusion that the spatial dimensions of a two-baryon bound state are to be 
large: R > 1 fm. The latter inequality translates into the limitation on the 
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bound-state momentum k ~ 1/R < 200 MeV. The correlation function C 
noticeably deviates from 1 for the values of q ~ k (see the figure). Strong 
correlation in the ir — n system are known to occur at the effective mass value 
~ 500 -7- 600 MeV, i.e. at the values of q in the same momentum range. One 
can conclude that in the q range of interest (q « 200 MeV) the correlations 
between pions are weak thus justifying the approximation made. 

One more implicit assumption made consists in neglecting the 7r-meson 
rescattering on the (heavy) particles. Such rescattering could generate some 
phase of the pion wave function ([H])- The magnitude of the rescattering 
effect can be roughly estimated as follows. Let a be a cross section of the 
emitted 7r-meson rescattering on the other baryon. The probability of the 
process is controlled by the ratio <j/(4ttR 2 ). In the considered momentum 
range \p\ = \q\/2 ~ 100 -j- 200 MeV the cross section of the nN scattering 
o^n k, 50 -T- 100 mbarn M. Then the rescattering parameter for R > 2 fm 
proves to be small < 25%. 

Summarizing, possible shallow bound states of two baryons can be inves- 
tigated by means of the pion interferometry. The momentum dependence of 
the two-particle correlation function in the "back-to-back" kinematics is con- 
trolled by the binding-energy value of two baryons. The latter can be deter- 
mined experimentally from the slope of the two-particle correlation function 
at (small) momenta of the 7r-mesons. The proposed method suits best for 
the search of possible (loosely) bound states of two resonances which binding, 
because of the large width, can hardly be observed as a shift of a peak in 
corresponding effective-mass distributions. 

The author is grateful to M.Polyakov and O.Patarakin for helpful discus- 
sions. 
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